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Abstract. We consider the response of a chaotic cavity in d dimensions to 
periodic driving. We are motivated by older studies of one-body dissipation 
in nuclei, and also by anticipated mesoscopic applications. For calculating the 
rate of energy absorption due to time-dependent deformation of the confining 
potential, we introduce an improved version of the wall formula. Our formulation 
takes into account that a special class of deformations causes no heating in the 
zero-frequency limit. We also derive a mesoscopic version of the Drude formula, 
and explain that it can be regarded as a special example of our calculations. 
Specifically we consider a quantum dot driven by electro-motive force which is 
induced by a time-dependent homogeneous magnetic field. 



1. Introduction 

The dynamics of a particle inside a cavity (billiard) in d = 2 or 3 dimensions is 
major theme in studies of classical and quantum chaos. Whereas the physics of time 
independent chaotic systems is extensively explored, less is known about the physics 
of time-dependent chaotic systems. The main exceptions are the studies of the kicked 
rotator and related systems jjj. However, the rotator (with no kicks) is a ID integrable 
system, whereas we are interested in chaotic (2D or 3D) cavities. 

Driven cavities were of special interest in the studies of one-body dissipation in 
nuclei fi 0, 0, |j. A renewed interest in this problem is anticipated in the field of 
mesoscopic physics. Quantum dots can be regarded as small 2D cavities whose shape 
is controlled by electrical gates. Another variation is driving a quantum dot by time- 
dependent magnetic field. In Section 6 we will explain that the calculation of the 
system response in the latter case can be regarded as a special example of the study 
in this paper. A similar observation applies to the case of a quantum dot driven by a 
homogeneous time-dependent electric field. However, in the latter case it is essential 
to take screening into account Q , and therefore our calculations no longer apply. 

We consider a system of non-interacting particles inside a cavity whose walls can 
be deformed. We define a single parameter x that controls this deformation. We 
would like to consider the case where x(t) — As'm(ilt) is being changed periodically 
in time, where A is the amplitude and £1 is the driving frequency. In particular we are 
interested in the small frequency limit, meaning fl <C 1/t co1 . Here r col is the typical 
time between collisions with the moving walls of the cavity. 

We will be interested in general deformations which need not preserve the billiard 
shape nor its volume. We can specify any deformation by a function -D(s), where s 
specifies the location of a wall element on the boundary (surface) of the cavity, and 
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D(s)Sx is the normal displacement of this wall element. There is a restricted class 
of deformations that are shape-preserving: they involve translations, rotations and 
dilations of the cavity. We will see that this class has special properties. Note that 
translations and rotations are also volume-preserving, in which case the associated 
time-dependent deformations can be described as 'shaking' the cavity. 

What is the rate in which the 'gas' inside the cavity is heated up? The answer 
depends on the shape of the cavity, the deformation D(s) involved, as well as on the 
amplitude A and the driving frequency f2. Also the number of particles N and their 
energy distribution p(E) should be specified. 

For non-interacting particles the solution of this problem is reduced to the analysis 
of one-particle physics. This observation is self-evident for non-interacting classical 



particles, but it is also true for non- interacting fermions (see Appendix A ). We would 
like to work within the framework of linear response theory (LRT) . In such case one 
can write 

±(H) = p(n) x ±(Anf (1) 

where the dissipation coefficient yu,(S7) is amplitude independent. The small-fi version 
of this formula can be written as 

cTt {n) = flV2 (2) 

where p, = fi(fl — > 0) is known as the friction coefficient, and V = Af2/v2 is the 
average root-mean-square (RMS) deformation velocity. For convenience let us define 
x as having units of length, such that V characterizes the velocity of the moving walls. 

A necessary classical condition for the validity of LRT is V -C v E where 
v E = (2E/m) 1 / 2 is the velocity of the particle |?], ||. We also assume that the motion 
of the particle inside the cavity is globally chaotic, meaning no mixed phase space M|. 
The criteria for having such a cavity are discussed i n |TT|, |l2] ] . The justification of LRT 
in the quantum-mechanical case is more subtle || |13|, and does not constitute a 
theme in this paper, although we do connect with the quantum case in Section^. The 
theory to be presented assumes that LRT is a valid formulation of the problem. 



As explained in Appendix A , LRT implies that the dissipation coefficient /i(f2) 
is related via a fluctuation-dissipation (FD) relation to a spectral function C E (f2). 
Namely, 

m = \j~p{E)dE J_ A [ ff(jB) c E (ft)] (3) 

Here p(E) is the energy distribution of the particles, and g(E) is the density of states. 
The function C e (oj) is the noise power spectrum of the generalized 'force' associated 
with the parameter x. This function is the main object of the present study, and its 
precise definition is in Section |^. We shall chiefly explore how C e (lj) depends on the 
type of deformation involved, but also discuss effects due to the cavity shape. 

In particular we are interested in the small frequency limit where p is related to 
the fluctuations intensity 

/oo 
C E (r)dr. (4) 
-oo 

The simplest estimate for i/ E , which we are going to call 'white noise approximation' 
(WNA) , leads (in case of a 3D cavity) to the well known 'wall formula' U 



Mb 



N f 

—mv E j> D(s) 2 ds (5) 
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where the subscript e implies that we are considering a microcanonical ensemble p(E), 
the number of particles is N, and the volume of the cavity is V. The above version 
of the wall formula has been derived for the purpose of calculating the so-called one- 
body dissipation rate in nuclei. The original derivation of this formula is based on a 
simplified kinetic picture Q. For alternate derivations using the LRT approach see 
H . For the generalization to any dimension d using the LRT-FD strategy see and 
further references therein. 

Our main purpose is to introduce an improved version of the wall formula, and 
to analyze the frequency dependence of /j.(S1) . This will involve a demonstration []l6f 
that for special types of deformations (namely dilations, translations and rotations) 
the small-O dissipation rate is remarkably different from the naive expectation. As an 
application, the mesoscopic version of Drude formula for the conductance of a quantum 
dot in a uniform time-dependent magnetic field reduces to the the calculation of C e (lu) 
for one of these special deformations (namely rotation), and leads to (see Section ||) 



where A is the area of the dot. 

For our improved wall formula, we show that it is essential to project out the 
special components of a general deformation, and only then to estimate v E using the 
WNA. If the assumption of strong chaos cannot be justified, further corrections are 
required due to correlations between successive bounces. 

The effect of interaction between the particles is not discussed in this paper. If 
the mean free path for inter-particle collisions is large compared with the size of the 
cavity then we expect that our analysis still applies. If the mean free path is much 
smaller, then we get into the hydrodynamic regime. In the latter case we have a drag 
effect, and the dissipation rate is determined by the viscosity of the gas via Stokes' 
law. 

2. The model system 

Consider a particle whose canonical coordinates are (r, p) moving inside a cavity. The 
Hamiltonian is 



where U(r) is the confining potential. We have introduced a (unitless) deformation 
'field' D(r), and x is the controlling parameter. In this paper we assume that U(r) = 
inside the cavity. The volume of the cavity is V. Outside the cavity the potential U(r) 
becomes very large. To be specific, one may assume that the walls exert a normal 
force /, and we take the hard wall limit / — > oo. With the above assumptions about 
U (r) it is clear that the deformation is completely specified by the boundary function 
D(s) = n(s) ■ D(s), where n(s) is an outwards unit normal vector at the boundary 
point s. 

Most of our numerical tests will refer to the 2D cavity illustrated in Fig. |l|a. 
It is a generalized two-dimensional Sinai billiard formed from concave arcs of circles 
with two different radii. Typical parameters used are a=2, 6=1, 01=0.2, 6*2=0.5, for 
which the average collision rate with the wall is (1/t w ) 0.63. This billiard has 
been chosen because it has 'hard chaos': There is no mixed phase space, and there 
are no marginally-stable orbits (see Section R). In Fig. [lib we show three example 




(6) 



H(r, p; x) = p 2 /2m + U(r - jcD(r)) 



(7) 
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Figure 1. a) The generalized two-dimensional Sinai billiard which has been 
used for our numerical studies, b) Three example deformations are illustrated. 
Note that they are shown exaggerated in strength. 



deformations. For illustration purposes we have selected three 'localized' deformations. 
See Tables 1 and 2 for a full list of deformations that have been tested in our numerical 
work. 

Associated with the parameter x is the fluctuating quantity 
&H(v, p; x) 



dx 



(8) 



where the time-dependence arises from that of the trajectory (r(i), p(t)). This quantity 
can be thought of as the generalized time-dependent 'force' associated with the 
parameter x. For the Hamiltonian (^) we can write 

T(t) = D(r) • VC/(r) = -D(r) • p. (9) 

Recognizing p as the force on the gas particle, we see that F{t) is a train of spikes 
(see Fig. ^]a). Namely, the fluctuating force T(t) consists of impulses whose maximum 
duration is tq = 2mv E / f. In the hard wall limit tq — > 0, and we can write 



' 2mv E cos(0j) Di S(t — ti) 



(10) 



where i labels collisions: ti is the time of a collision, Di stands for D(si) at the 
location Sj of a collision, and v E cos(9i) is the normal component of the particle's 
collision velocity. The above sequence of impulses is characterized by an average rate 
of collisions l/r col . The quantitative definition of r col is postpone to Section |[ Note 
however that t co1 may be much larger that the ballistic time r bl . The ballistic time is 
the average time between collisions with the boundary. We have t co1 3> r hl whenever 
a deformation involves only a small piece of the boundary. Finally we note that if 
the deformation is volume-preserving then {T(t)) = 0. Otherwise it is convenient to 
subtract the (constant) average value F(x) from the above definition of T{t). This 
convention is reflected in our illustration (Fig.^a). 

We define the auto-correlation function of T(t) as follows: 

C E (r) = {T{t)T{t + r)) E (11) 

The subscript e, whenever used, suggests that the average over initial conditions is 
of microcanonical type, with energy E. Note that C e (t) is defined using the time 
independent ('frozen') Hamiltonian, and therefore is independent of t. The auto- 
correlation function C E (r) can be handled as a time-average rather than an ensemble- 
average (by ergodicity). The resulting construction is illustrated in Fig. ^b, where we 
illustrate the projection of T{ti)T(t2) onto the t = £2 — t\ axis. The contribution 
for the self correlation is shaded. The forms of the resultant C E (r) and its Fourier 
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transform C E {uo) = J C E (r) exp(iuiT)d,T are illustrated schematically in Figs. |^c and 
||d. Note that the u> — > limit of C^w) is equal to the area under C(r). 

The auto-correlation function C E (r) consists of a r = ('self') peak due 
to the self-correlation of the spikes, and of an additional smooth ('non-self') 
component due to correlations between successive bounces. This implies JsJ that 
pronounced correlations are usually characterized by the time scale r bl , rather than r col . 
Consequently the associated frequency scale for non- universal structures is to ~ l/r bl . 
Another relevant time scale is the ergodic time i org which is the inverse of the average 
Lyapunov (instability) exponent. Beyond i org the correlations become vanishingly 
small. Non-negligible tails may arise only if the motion has marginally stable orbits. 

As explained in the Introduction, we shall be most interested in the noise intensity 
v E defined by (|4|). Observing that T(t) is linear in D(s), it follows that the noise 
intensity is a quadratic functional 

v E = //efei<2s2-D(si)7 E (si,s 2 ).D(s 2 ), (12) 

where the kernel 7 E depends on both the cavity shape and the particle energy E ||. 
Furthermore, billiards are scaling systems in the sense that a change in E leaves the 
trajectories unchanged. From this and ( |To| ) we have the scaling relation 7 E (si,S2) = 
m 2 v E ■ 7(si, S2), where the scaled kernel depends entirely on the geometrical shape of 
the cavity. However, the reason for being interested in approximations for v E is that 
the exact result for the kernel 7 is very complicated to evaluate, and involves a sum 
over all classical paths from Si to s 2 (see [||). 

3. Quantum-classical correspondence 

This paper applies classical physics in order to analyze the response of a wide class 
of systems, including mesoscopic systems where quantum mechanics may play a role. 
How much of a compromise is a classical analysis of the dissipation? This question 
has been addressed in (r| ||. In the level of one-particle physics the answer is 
as follows: within the framework of LRT the only difference between the classical 
formulation and the quantal one is involved in replacing the classical definition of 
C e (t) by the corresponding quantum- mechanical definition. In the level of many (non- 
interacting) particles the only further modification is associated with the application 



of the FD relation, as discussed in Appendix A (See Eq.( A. 6 )). We would like to 
re-emphasize that we assume in this paper that we are in an (A, f2) regime where LRT 
is a valid formulation. The quantum adiabatic regime (extremely small f2), and the 
non-perturbative regime (see discussion in Jl4f| ) are excluded from our considerations. 

Thus the only remaining question is whether a classical calculation of C E (ui) 
is a good approximation quantum-mcchanically. The answer is that the quantum- 
classical correspondence here is remarkable. It has been tested for a few example 
systems §]], In Fig. | we demonstrate correspondence for the stadium billiard 

for three types of deformations: DI (dilation), W2 (periodic oscillation around the 
perimeter), and P (wide 'piston' existing only on the top edge). The RMS estimation 
error is 3% for the classical calculation and 10% for the quantum calculation. The 
quantum estimate of C E (w) amounts to computing boundary overlap integrals of the 
eigenfunctions (see Jl6|]). We have used all 451 states lying between wavenumbers 
398 < k < 402, where the mean level spacing is A w 8.8 x 10~ 3 in co units. Note 
that there are ~ 10 2 de Broglie wavelengths across the system. The stadium was 
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Figure 2. The fluctuating force looks like a train of impulses (a). Due to 

ergodicity the autocorrelation function C(t) can be regarded as a time average (b). 
The resultant autocorrelation function (c) and the associated power spectrum (d) 
may be characterized by non universal features. See text for further explanations. 



chosen because it enables efficient quantization using the method of Vergini and 
Saraceno pol. An especially good basis set is known for this shape H. 



4. The white noise approximation 

The most naive estimate of the fluctuations intensity is based on the WNA. Namely, 
one assumes that the correlation between bounces can be neglected. This corresponds 
H to the local part of the kernel (12). In such case only the self-correlation of the 
spikes is taken into consideration and one obtains j8| 

is E « (2mu E ) 2 ^cos 2 (^) D\ 6(t-ti)J (13) 

and from here (see ^]) using ergodicity, 

u E « 2m 2 vl(\cos8\ 3 )^ <£[D(s)] 2 ds (14) 

where the geometric factor for d = 2, 3, • • • is (| cos(#)| 3 ) = 4/(37r), 1/4, • • •. If we can 
use the convention |D(s)| ~ 1 over the deformed region (and zero otherwise), then 
we can write the WNA as u E = (2mv E ) 2 x (l/r col ) where (1/t co1 ) defines the effective 
collision rate. For a more careful discussion see Appendix F of Ref.f|. Note again 
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Figure 3. Agreement between quantum and classical Ce(o;) in the two- 
dimensional quarter-stadium billiard for three example deformations (see text). 
In each case classical is shown as a thick line, and quantum a thin line. The 
j/-axis has been displaced to clearly show the w — » behavior. The singular peak 
at ui = 7r is due to the 'bouncing ball' orbit. 



that r col can be much larger than the ballistic time r bl in the case that only a small 
piece of the boundary is being deformed. 

The use of the WNA can be justified whenever successive collisions are effectively 
uncorrelated. The applicability of such an assumption depends on the shape of 
the cavity (which will determine the decay of correlations via the typical Lyapunov 
exponent) as well as on the type of deformation involved. If we have the cavity 
of Fig. Ilk, and the deformation involves only a small piece of the boundary (e.g. 
see Fig . Qjb) , then successive collisions with the deformed part of the boundary are 
effectively uncorrelated. This is so because there are many collisions with static pieces 
of the boundary before the next effective collision (with non zero Di) takes place. 
If the deformation involves a large piece (or all) of the boundary, we can still argue 
that successive collisions are effectively uncorrelated provided D(s) is 'oscillatory' 
enough (ie changes sign many times along the boundary). These expectations are 
qualitatively confirmed by the numerical results of Fig. |5| Here we show a sequence 
of deformation types for which the WNA performs increasingly well: FR (for which 
sensitivity to the vertical least-unstable periodic orbit causes large correlation effects 
and large deviations from WNA), W8 (oscillatory deformation changes sign many 
times around the perimeter, giving better agreement with WNA), PI (localized 
'piston' type deformation, for which WNA is good), and DF (random function of 
zero correlation-length along the perimeter, showing complete WNA agreement). 

The numerical evaluation of C e (lo) throughout this work is performed by squaring 
the Fourier transform of a single long sample of T(t) ( ~ I0 6 consecutive collisions). 
Ergodicity ensures that the properties of a single trajectory reproduce the desired 
ensemble-average (• • -) E . In practice the power spectrum of a single sample is a 
stochastic quantity with no correlations in w-space. To estimate the underlying noise 
spectrum C e (uj) a smoothing convolution in cj-space is performed. In the figures a 
smoothing width of I0~ 2 is typical, giving 3% RMS estimation error. The <5-function 
nature of T(t) is handled by convolving in the time-domain with a suitably-narrow 
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Figure 4. The white-noise approximation estimate (WNA is the horizontal 
dotted line) compared to actual C E (t<j) power spectra for some example 
deformations of the 2D generalized Sinai billiard, with m = v = 1. (The RMS 
estimation error of 3% can be seen as multiplicative noise with short correlation 
length in u>). Deformation functions are defined in Table I, and discussed further 
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key 


description 


surface deformation function D(s) 


CO 


constant 


1 


Wn 


n periods 


cos(27ms/L) 


DF 


diffuse 


random[-l,l] (equivalent to Woo) 


FR 


fracture 


sgn(x(s)) if on top or bottom, else 


SX 


shift-x 


sgn(x(s)) if on left or right, else 


PI 


piston 1 


10exp(-|a 2 ), a=(s/L- 0.3) /0.01 


P2 


piston 2 


10exp(-±a 2 ), a=(s/L - 0.6)/0.005 


WG 


wiggle 


5aexp(-±a 2 ), a=(s/L - 0.25)/0.02 



Table 1. Key to deformation types used for numerical 2D billiard experiments 
in this paper. L is the billiard perimeter. The deformation is described by a 
function D(s), where s is measured counter-clockwise along the perimeter with 
s = at the upper left corner. In the 'fracture' and 'shift-x' cases we use the 
horizontal Cartesian coordinate x(s). 

Gaussian. This enables the signal to be sampled uniformly in time, and hence we can 
benefit from use of the Fast Fourier Transform procedure. 

It might be asked whether the exponential growth in sensitivity to numerical 
round-off error invalidates the computation of the properties of a long classical 
trajectory. The answer is no: it has been shown that in simple two-dimensional 
chaotic maps such as ours, a numerically-generated 'pseudo-trajectory' shadows (is 
very close to) a true trajectory with slightly different initial conditions [^]. However, 
as we shall see, the differences in to — > behavior (in the hard chaos case) do not in 
fact rely on correlation properties over times any longer than i org . 
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Figure 5. The WNA estimate compared to actual C(oj) noise power spectra 
for example 'special' deformation types: DI (dilation), TX (translation) and RO 
(rotation). See Table II for definitions. The WNA fails to predict the vanishing 
in the small u> limit. 

key description deformation field 

DI dilation about origin D(r) = r 

TX x-translation D(r) = e x 

TY ^-translation D(r) = e y 

RO rotation about origin D(r) — e z x r 

Table 2. Key to the four 'special' deformations in 2D. The unit vectors e x and 
ej, are in the plane (see Fig. 1), and e z is in the perpendicular direction. In the 
case of dilation and rotation D could be made unitless by dividing by a constant 
length. 




5. 'Special' deformations 

The WNA dramatically fails (see Fig. ||) for dilation, translations and rotations (see 
Table |^ for their definitions in 2D). It is not surprising that the WNA is 'bad' for 
these deformations because their D(s) are slowly-changing delocalized functions of s. 
However, what is remarkable is that C E (u>) for this type of deformations vanishes in the 
limit lu — ► 0. Such deformations we would like to call 'special' More generally, we 
would like to say that a deformation is 'special' if the associated fluctuation intensity 
is u E = 0. 



A result that follows from the considerations of Appendix B is that a linear 
combination of special deformation is also special. Therefore the special deformations 
constitute a linear space of functions. We believe that this linear space is spanned by 
the following basis functions: one dilation, d translations, and d(d — l)/2 rotations. 
However we are not able to give a rigorous mathematical argument that excludes 
the possibility of having a larger linear space. In other words, we believe that any 
special deformation can be written as a linear combination of dilation, translations 
and rotations. 

We will explain the observed v E = 0, starting with the case of translations and 
dilations. For translations we have D = e, where e is a constant vector that defines a 
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direction in space. We can write J-(t) = (d/dt) 2 G(t) where Q(t) = —me ■ r. A similar 
relation holds for dilation D = r withC/(£) = — \mv 2 . It follows that C{uj) — u> 4 Cg(u), 
where Cg(w) is the power spectrum of Q{t). If Cg(w) is a bounded function (as it 
must be when correlations are short-range), it immediately follows that C(0) = 0. 
Moreover since Q(t) is a simple function of the particle position, we can assume it is 
a fluctuating quantity that looks like white noise on timescales > t or(5 . It follows that 
C(u>) is generically characterized by uj behavior for cither translations or dilations. 

We now turn to consider the case of rotations. This case is of particular interest 
because of its relation to the Drude conductance calculation in a uniform driving 
magnetic held (see the following section). For rotations we have D = e x r, and we 
can write !F(t) — (d/dt)Q(t), where Q{t) = — e ■ (r x p), is a projection of the particle's 



angular momentum vector |15. Consequently C(uj) = uj 2 Cg(uj)- Assuming the 
angular momentum is a fluctuating quantity that looks like white noise on timescales 
> t crs , it follows that C(oj) is generically characterized by uj 2 behavior. 

Thus we have predictions for the power-laws in the regime u> < 1 / t OIg for special 
deformations (assuming hard chaos). These have been verified numerically in our 
previous work |1^| , with a special emphasis on the case of dilation. The case of dilation 
plays a vital role in a highly-successful numerical billiard diagonalization method that 
has been introduced recently [[l9) . 

For special deformations we have C(uj) — in the limit uj — 0, and consequently 
the dissipation coefficient vanishes (/i = 0). It should be noted that for the case 
of a general combination of translations and rotations this result follows from a 
simpler argument. Taking fl — > while keeping Ail constant corresponds to constant 
deformation velocity (i =const). Transforming the time-dependent Hamiltonian into 
the reference frame of the cavity (which is uniformly translating or rotating with 
constant velocity) gives a time-independent Hamiltonian. In the new reference frame 
the energy is a constant of the motion, which implies that the system cannot absorb 
energy (no dissipation effect), and hence we must indeed have /i = 0. 

6. Drude mesoscopic conductance for 2D dot 

Consider a 2D quantum dot in a homogeneous (perpendicular) magnetic held (see 
Fig. ||b). The one-particle Hamiltonian is 

W(r,p;$(t)) - -^[p-eA(r;$(t))] 2 + C/(r) (15) 

The dot is defined by the confining potential U(r), and we choose the magnetic field 
as the controlling (driving) parameter. Periodic driving means $(i) = Asin(ftt). The 
vector potential is given by 

A(r;#) = i(?a) xr (16) 



2 VA 

where A is the area of the dot, $/A is the magnetic field, and z is its (perpendicular) 
direction. 

Referring to Eq.(2) one should realize that by Faraday's law V = $ is the 
induced electromotive force (measured in volts). Hence fi is just the conductance. 
The fluctuating quantity that is associated with $ has the meaning of electric current: 

X(t) = = — (z x r) ■ v (17) 

W 9$ 2A V ' y ' 
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In the conventional ring geometry (Fig. ^a) the current is just I{t) = (e/L)v, where L 
is the perimeter, and v is the tangential velocity. In the general cavity case (Fig. |6|b) 
I(t) can be thought of as the angular momentum of the charge. 

The Drude mesoscopic conductance is given by the frequency-dependent version 



of the FD relation Eq.(A.E). With one-particle density of states corresponding to 2D 
gas Eq.( [A.6| ) becomes 

/x(O) = -^<7z(n) (18) 

where the Fermi velocity is related to the Fermi energy E F — \mv\. The power 
spectrum of the electric current Cx(ui) is the Fourier transform of the current-current 
correlation function. In standard derivations of the Drude formula it is assumed that 
this correlation function is exponential: 

Cz(r) ~ ^vl exp (19) 

leading to the Lorentzian Eq.(||). However, for a given dot shape Cj(t) is not really 
an exponential, but rather reflects the system-specific geometry. Below we discuss two 
limits in which we can obtain approximations for Cx(r) and hence (via Eq.(|lS|)) for 
the frequency-dependent conductance /x(f2). 

The current 1(t) is a piecewise constant function of time. It is constant between 
collisions with the walls because of conservation of angular momentum. The derivative 
of this quantity, !F(t) = I, is a train of spikes. It formally coincides (using @) with 
the 3-{t) of the deformation D(r) = (e/(2mA))z x r, corresponding to rotation around 
the z axis. It follows that the current-current correlation C%(t) is trivially related to 
the T{t) correlation function C(r) as follows: 

<3r(w) = \c{u) (20) 

Thus we see that the calculation of 'conductance' is formally equivalent to a special 
case of deformation, namely a rotation. 

There are two limits in which we can get approximation for Cj(lu). For small 
frequencies u> <C (1/t ooJ ) we may use the following simple estimate: 

CxH = <X 2 )x2r COI = I e 2 ^ w 2 xr col (21) 

The first equality can be taken as an operative definition of the correlation time r col in 
the context of this calculation. Obviously, up to a system-specific geometrical factor 
this result (~ uj°) agrees with the standard Drude result. For ring geometry one should 
make the replacements (r 2 ) i— > (L/(27r)) 2 and A i— > 7r(L/(27r)) 2 where L is the length of 
the wire (perimeter of the ring). Thus one obtains Cx(w) = ((e/L)u F ) 2 leading to the 
standard- looking Drude formula for a mesoscopic wire \i = (iV/L 2 ) x (e 2 /m) x r col . 

In the limit u> 3> (1/t co1 ) we can get a much more satisfying result. The fluctuating 
quantity T{€) = i(t) is the same as (0) with D(r) = (e/(2mA))z x r, corresponding 
to rotation. Using the WNA of Eq. (CT) , and dividing by oj 2 as in j20|) we get 



Cx(w) 



2 



i 

"^2 



(22) 



Again, up to system-specific geometrical factor this result (~ u> 2 ) agrees with the 
standard Drude result. The latter expression should become exact as we go to large 
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Figure 6. Two possible mesoscopic geometries which exhibit conductance when 
driven by a magnetic field: a) conventional ring of perimeter L enclosing the time- 
dependent flux, b) ballistic two-dimensional chaotic dot (cavity) of area A in a 
uniform time-dependent magnetic field. 



frequencies, where the only significant contribution comes form the self-correlation of 
the T(t) spikes (see Fig. 0d). 

Eq.(pl|) leads (via (18)) to the small frequency Drude result, while the WNA of 
Eq.(^2|) gives the Lorentzian tail of the Drude result. An exact result for the frequency- 
dependent conductance can be calculated numerically for a given geometrical shape. 
In Fig. [?] we display a plot of /x(f2) oc Cx(fi), which shows both the constant behavior 
at small fi and the convergence to the large-f2 WNA approximation. System-specific 
features are expressed by the deviation from standard Lorentzian in the intermediate 
frequency regime. 

Finally we consider driving a quantum dot with homogeneous electric field in the 
x direction, in which case the Hamiltonian contains the interaction term —e£(t)x. 
For calculation of the response in such a case one should evaluate the dipole-dipole 
correlation function C-p(r) where V(t) = ex. The latter is related to translations, 
where the deformation field is D = x. Consequently we get C-p{r) = (1/cj 4 )C(t). 
However, this result is not of great interest, because the screening effect leads to 
modification of the effective one-particle Hamiltonian, such that the actual electric 
field inside a quantum dot is much smaller than the applied field. 



7. The white noise assumption revisited 

In Section [| we have assumed that generic fluctuating quantities such as r 2 and e • r 
and e • (r x p) have a white noise power spectrum for u -C 1/t w . In section ^ we 
are going to suggest that this white noise assumption is approximately true for any 
fluctuating quantity T(f) that comes from a normal deformation (the term 'normal' 
will be defined there). 

Obviously, the goodness of the 'white noise assumption' in the two cases 
mentioned is related to the chaoticity of the system, and should be tested for particular 
examples. This has been done so for the cavity of Fig. [I] (see p6|, and Figs. [I] and ^|). 
This cavity is an example of a 'scattering billiard' and so exhibits strong chaos O]. If 
the motion is not strongly chaotic we may get a C(r) that decays like a power law (say 



1/r 1 7 with < 7 < 1) rather than an exponential O, |23|, 24, g5| £6|]. In such case 
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Figure 7. Calculation of dissipation coefficient fJ.(Q) (arbitrary units) for driving 
of a chaotic mesoscopic billiard system with a constant magnetic field at frequency 
S7. The billiard chosen is the 2D generalized Sinai of Fig. lla. The dotted line 
WNA(JF) is the high-frequency estimate assuming JF(t) = 2is white noise. The 
convergence to this oj -2 result is clearly visible in the log-log inset plot. 




the universal behavior is modified: we get o-> -7 behavior for C e (lo) at small frequencies 
(v E diverges), signifying faster-than-diffusive energy spreading in Eq^AJ!) |2(|. The 
stadium is an example where such a complication may arise: an ergodic trajectory 
can remain in the marginally-stable 'bouncing ball' orbit (between the top and bottom 
edges) for long times, with a probability scaling as t^ 1 |^3|, ^4, E5 . Depending on the 
choice of D(s) this may manifest itself in C(r). For example, in Fig. || the deformation 
P involves a distortion confined to the upper edge, and the resulting sensitivity to the 
bouncing ball orbit leads to large enhancement of the fluctuations intensity C(w=0), 
and is suggestive of singular behavior for small lu. 

If the billiard has a mixed phase space (which is the generic case), then the 
integrable component does not contribute to diffusive energy spreading. Proposals 
have been made to account for this via a phase-space volume factor |l(| . 



8. Decomposition of general deformations 

The failure of the WNA for 'special' deformations also extends to the much wider 
class of deformations which are similar to special. This is demonstrated in Fig. |[ It 
should be emphasized that this failure happens even if the cavity is strongly chaotic. 

We seek an analytical estimate for C(w), and in particular for its zero- frequency 
limit v. This estimate should apply to any (general) deformation, including the case of 
'close-to-special' deformations. It would be useful to regard any general deformation 
as a combination of 'special' component and 'normal' component. The formulation 
of this idea is the theme of the present section. Supporting numerical evidence is 
gathered in the next section. 

The special deformations (for which we have v = 0) constitute a linear space, 
meaning that any sum of special deformations is also a special one. Now we would 
like to conjecture that there is also a linear space of 'normal' deformations. By 
definition, for 'normal' deformation T(t) looks like an uncorrelated random sequence 
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Figure 8. The failure of the WNA estimate for C(uj) for deformation types CO 
(similar to DI) and SX (similar to TX). The WNA is clearly a vast overestimate 
of the small-ij limit. See Tables I and II for explanation of deformation types. 



of impulses, and consequently the WNA is a reasonable approximation. The notion of 



randomness can be better formulated as in Appendix C leading to Eq.(C4). However 
in practice (|C.4[ ) is not useful, because it cannot be applied as an actual classification 
tool. (Eq.(CT) is never satisfied exactly). Still we are going to demonstrate that 
there is a unique way to identify the subspace of normal deformations, if we insist on 
a maximal (i.e. the most inclusive) definition of this subspace. 

It is important to clarify the heuristic reasoning of having linear space of normal 
deformations. The Tit) that corresponds to some normal deformation D(s) looks 
like white noise. It means that only self-correlations of its spikes are statistically 
significant. If we have two such generic quantities, say Ti(t) and T 2 (t), then we 
expect T\ (t) + T 2 (t) to share the same property. 

The correlation function of T(t) = T\ (t) + T 2 (t) can be written formally as 

Ci +2 (r) = Ci(r) + C 2 (r) + 2C 1)2 (r) (23) 
where C\^ 2 {t) is the cross-correlation function. In |Appcndix B we argue the following 

C li2 (r)dT = 

if l=general, 2=special (24) 



This result is exact, and does not involve any approximation. In Appendix C we argue 
the following 



Ci, 2 (t) ps c x 
if l=normal, 2 



<j^D 1 {s)D 2 (s)ds 



5(r) 



=general (25) 

where c = 2m 2 v^(\ cos6*| 3 )/V. This result is an approximation, which is expected to 
be as good as our assumption regarding the 'normality' of the deformation Di(s). 
Consider now the case where -Di(s) is normal and D 2 (s) is special. Both Eq.(p4|) 
and Eq.d2a) should apply. But these equations are consistent if and only if Di(s) is 
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orthogonal to D 2 (s). We say that D 1 (s) and D 2 (s) are orthogonal (1 _L 2) using the 
following definition: 

orthogonality <j> D 1 (s)D 2 (s)ds = (26) 

Thus we have proved that normal deformations must be orthogonal (in the sense of 
(p6|)) to special deformations. Obviously we have proved here a necessary rather than 
a sufficient condition for 'normality'. However, if we insist on a maximal definition 
for the subspace of normal deformations, then we get a unique identification. Namely, 
a deformation is classified as 'normal' if it is orthogonal to the subspace of special 
deformations. 



The practical consequences of Eq.(24) and Eq.(E5|) are as follows: 



^1+2 = v\ if l=general, 2=special (27) 

and 

^1+2 » V\ + v% + 2c j> D 1 (s)D 2 (s)ds 

if l=normal, 2=gcncral (28) 
These results are tested in the next section. 



9. Addition of deformations: numerical tests 

On the basis of the discussion in the previous section we define normal deformation 
as those that are orthogonal to all special deformations, in the sense of Eq. (E6|). 
Obviously there are 'good' normal deformations for which the WNA is an excellent 
approximation (PI and W8 in Fig. [|, for example), and there are 'bad' normal 
deformations for which the WNA is not a very good approximation (FR in Fig. [|, 
and the normal component in Fig. |l4|b). In this section we present numerical evidence 
that verifies the theoretical results of the previous section, and investigate how 'bad' 
a normal deformation has to be for them to break down. 

From what we have claimed it follows that if Di(s) and D 2 (s) are orthogonal 
normal deformations, then ^1+2 = ^1 + v 2 . We could as well write 

C 1+2 (lo) » Ci(w) + C 2 (w) 

if l=normal, 2=normal, and 1 _L 2 (29) 

because by assumption the three correlation functions are approximately flat. We 
demonstrate this addition rule in the case of two 'good' deformations which are 
orthogonal in Fig. ^. We found that small 'pistons' (P2 is significant on only ~ 1/50 
of the perimeter) were needed to achieve addition of the accuracy (a few %) shown. 
However, the restriction on the 'wiggle' type of deformation was somewhat more 
lenient (WG is ~ 5 times wider than P2 yet obeys the WNA better than P2 does). 

In general we observe that the quality of the addition rule is limited by the 
deviation from the WNA of the better of the two deformations. In Fig. [Io| we see 
that if both -Di(s) and D 2 (s) are bad, then also the addition rule ( p9f ) becomes quite 
bad. Fig. [ll] shows that the addition rule (|2^) is reasonably well satisfied also if either 
D 1 (s) or D 2 (s) is a 'good' normal deformation. We have chosen -Di(s) as WG (good), 
and D 2 (s) as SX which is almost completely dominated by the special x-translation 
deformation. The addition rule (p9|) is obeyed at all to. This proves that our assertions 
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Figure 9. Addition of two 'good" normal deformations (1=P2, 2=WG). The 
two are orthogonal in the sense of (Bfj). That they are 'good' can be seen by their 
good agreement with their WNA results (horizontal arrows) . The power spectrum 
of the sum agrees well with the sum of the power spectra. 

Eq.(p5|) about the vanishing of Cf°2 BOlf ( r ) is indeed correct. It holds here as a non- 
trivial statement (^(s) is general and 'bad'). 

Finally, we consider the case where D\(s) is general and £*2(s) is special. This is 
illustrated in Fig. |lj. The addition rule ( p9| ) becomes exact in the limit of small 
frequency corresponding to the vanishing of Ci^{u — > 0) as implied by Eq. (Bjj) . 
In particular this implies that ^1+2 — v \- Note that there is no condition on the 
orthogonality of -Di(s) and Z?2(s)). This will be the key to for improving over the 
WNA, which we are going to discuss in the next section. 

In drawing the above conclusions it is important to note that symmetry effects 
can play a deceptive role if the cavity shape has symmetry (our example Fig. nlis 
in the Civ symmetry group) . In Fig. |l3| we demonstrate that the addition rule ( [29|) 
is very accurately satisfied at all uj if -Di(s) and £>2(s) belong to different symmetry 
classes of the cavity. Orthogonality of -Di(s) and D2(s) is not sufficient to explain 
this perfect linearity of addition of C E (u>). Rather, it follows from the symmetry 
of the kernel 7 E (s 1 ,s 2 ) of Eq.(|l|). The cross-terms in ( |l2[ ) rigorously vanish when 
such deformations are added. The consequence is that in order to demonstrate the 
assertions of this and of the previous section, we had to choose deformations of the 
same symmetry class, or which break all symmetries of the cavity. 

10. Beyond the WNA 

It is possible now to consider the case of general deformation, and to go beyond the 
WNA. Given a general deformation D(s) we should project out (subtract) all the 
special components, leaving the normal component, and only then apply the WNA. 
In Fig. |lj we demonstrate this decomposition for the deformation (CO + W16) and 
the deformation SX. 

The special deformations constitute a linear space which is spanned by the basis 
functions: one dilation, d translations, and d{d — l)/2 rotations. (For d=2 they are 
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Figure 10. Addition of two 'bad' normal deformations (1=FR, 2=SX). The 
two are orthogonal in the sense of That they are 'bad' is shown by a 

lack of agreement with their WNAs. The power spectrum of the sum is badly 
approximated by the sum of the power spectra (non-linear addition). 
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Figure 11. Addition of a 'good' normal deformation fl=WG) to a general 
deformation (2=SX). The two are orthogonal in the sense of (pd). The power 
spectrum of the sum agrees well with the sum of the power spectra. 



listed in Table II) . For a general cavity shape these basis functions are not orthogonal. 
However, because they are linearly independent, we can use standard linear algebra 
to build an orthonormal basis {-Dj(s)} of special deformations. The special (||) and 
the normal (_L) components of any given deformation D(s) are therefore 

#|| (s) =EiOiA(8) 

D ± (s) = D(s)-Z»||(s) (30) 
where the coefficients are 



a, = <j) D(s)Di(s) ds. 



(31) 
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Figure 12. Addition of a general deformation (1=FR) to 'special' deformation 
(2=TX). The power spectrum of the sum coincides with the sum of the power 
spectra in the limit uj — > 0, as implied by Eq. (B7f). 



C(co) 




Figure 13. Addition of two 'bad' general deformations which come from 
different symmetry classes of the cavity (1=W2, 2=FR). The two must also be 
orthogonal, by symmetry. The deviation from linear addition (solid line varying 
about zero) vanishes at all uj. 



The improved approximation for v applies the WNA only to the normal component, 
giving 



2m 2 «i(|cos0| 3 )^ Ids[D ± (s)f 



(32) 



which we name the IFIF (Improved Fluctuations Intensity Formula). In the particular 
case of d=3, substitution of this result into the microcanonical FD relation gives an 
'improved wall formula' consisting of the replacement of D(s) by D±(s) in Eq. (||). 

In Fig. [U we use the IFIF to estimate v for two examples. The first is a 
deformation (CO + W16) whose normal component is 'good', due its oscillatory 
nature. The deviation from a flat white power spectrum is ~ 20% for the normal 
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Figure 14. Decomposition of general deformations D(s) into orthogonal 
'normal' and 'special' components. The general deformation is CO + W16 in 
subfigure (a), and SX in subfigure (b). The naive_W NA Eq.(0) is indicated by 
short solid line. The improved (IFIF) result Eq.(|32[) is indicated by long dashed 
arrow. 



component. The IFIF result Eq.(|32|) is accurate to a few percent. It is a much better 
estimate of the actual v compared with the naive WNA Eq.(|l4|) which overestimates 
the correct value by a factor of 2.2. In the second example the deformation is SX. 
The resulting normal component is 'bad'. Its power spectrum fluctuates by a factor of 
about 10 in the lu range shown. Consequently the IFIF is limited in its accuracy, and 
the correct value for v is underestimated by a factor of 2.5. However, it is still a great 
improvement over the naive result Eq.(|14|). In this second example we can extract 
another prediction about C E (u). The special component is a factor ~ 10 larger than 
the normal component. Therefore the (J 2 behavior at small ui is almost entirely due to 
the 'rotation' component. The prefactor of the w 2 behavior need only be found once 
for each billiard shape (see Section ^). This saves computation and gives an extra 
information about the dissipation rate at finite driving frequency. 

A few concluding remarks regarding the history of the wall formula are in order. 
It has been known since its inception that the naive wall formula gives unphysical 
answers in the case of constant-velocity translations and rotations. This was first 
regarded as a kinetic gas 'drift' effect |2j. It should be noted that the recipe presented 
in 1 3], namely to subtract this drift component, is equivalent in practice to the recipe 
([30) that we have presented here, provided we ignore dilations. It is also important 
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Figure 15. a) A deformation of the stadium which moves the 'center of mass' 
(centroid) of the cavity to the right (from th e Ho t to the crosshairs symbol). This 
deformation is orthogonal (in the sense of Q26|) ) to all special deformations, in 
particular, all translations, b) An example volume-preserving deformation of an 
elongated approximately-rectangular cavity (/3 -C 1) which nevertheless has a 
large overlap with dilation. It can be shown that this results in an IFIF estimate 
of « 4/3 times that of the naive WNA. In both diagrams the undeformed shape is 
shown as a heavy line, the deformed one as a thin line. 
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to realize that the argumentation in |2j for this subtraction appears to be ad hoc, 
being based on a 'least-structured drift pattern' reasoning. A stated condition on this 
subtraction was that the resulting deformation preserve the location of the 'center 
of mass' (centroid) of the cavity, for reasons particular to the nuclear application 

This condition seems to have become standard practice in numerical tests of 
the wall formula [[l0[ ^7|, [2^]. However, as the system in Fig. |l5| a illustrates, 
this condition is generally not equivalent to the above subtraction of translation and 
rotation components p0[ |. This seems to invalidate the theorem presented in Section 
7.1 of S. Where the flaw in their reasoning lies we are not sure. 

The consideration of the special nature of dilations is absent from the literature. 
Even if we restrict ourselves to volume-preserving deformations (the case for the 
nuclear application), then deformations of certain cavities can be found for which 
the dilation correction is significant; we illustrate this in Fig. p"5|b. This correction can 
only be large if the cavity has a large variation in radius (i.e. is highly non-spherical). 
We suggest this as a possible reason why major discrepancies due to dilation have 
not emerged in the numerical tests of the wall formula until now. Such tests have 
generally been of shapes close to a 3D sphere ||, [2?], [2^, |29) . 

Hence we believe that the recipe we have presented, along with the associated 
theory and in conjunction with the particular power-law dependences, is a significant 
step in the treatment of one-body dissipation. 
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Appendix A. Linear response theory of dissipation 

Given a parametric Hamiltonian H(Q,P;x(t)), and given initial conditions, one 
defines the energy £(t) = H(Q(t), P(t); x(t)) and the fluctuating quantity T{t) = 
—dH/dx(Q(t), P(t); x(t)). With no approximation we have 

£ (t) - £(0) = f T{t')x{t')dt'. (A.l) 



Jo 

Using the same steps as in || one obtains the following result for the variance of the 
energy spreading: 

8E(tf = f f C E (t" ~t')F(t" -t')dt'dt" (A.2) 
Jo Jo 



where C E (r) is defined by Eq. (Jl l|) . Microcanonical averaging has been taken over 
the initial conditions. The function F(t) = (x(t)x(t + r)) is the velocity-velocity 
correlation of the driving. For periodic driving x(t) = Asm(£lt + phase) it is formally 
convenient to average over the initial phase and one obtains F(t) — \{AH) 2 cos(17r). 

For a chaotic system C e (t) is characterized by some correlation time r cl . For 
t ^> t c1 one obtains diffusive spreading 5E(t) 2 = 2D E t where the diffusion rate is 

D E = ^C E (Q) x \{Anf (A.3) 

which for small frequencies goes to D E = \v E V 2 , where v E = C(0) as defined in the 
Introduction. The picture to keep in mind is that of the fluctuating F(t) causing a 
random walk in energy space via Eq.(A.l) for times t ^> r cl . As explained in ^[ || 



the resulting diffusion in energy space implies systematic growth of the average energy. 
It is important to realise that this growth happens even if the random walk is locally 
unbiased: such is the case when changing the parameter x preserves the volume of a 
given energy-shell in phase-space. (For a deforming billiard system this corresponds 
to preservation of the billiard volume). The rate of energy growth is related to the 
diffusion as follows: 

where p(E) is the energy distribution of the particles, and g{E) is the one-particle 
density of states. The growth is therefore an effect of the ^-dependence of both the 
diffusion rate and the density of states. 

The rate of dissipation can be written as in Eq.([T| ) or as d{TL)/dt = [iV 2 in 



the small frequency limit. Combining this with Eq.(A.4) implies a relation between 
the dissipation coefficient /j, and the function C e (lu). The most familiar version of 
this FD relation is obtained for small frequency under the assumption of a canonical 
distribution p(E) oc g(E) exp(— E/^bT)), leading to 

" = »hf u (A - 5) 

where v should be calculated for a canonical distribution. This result should be 
multiplied by the number of non-interacting classical particles. 



The use of Eq.(A.4) can be justified also for non- interacting fermions [BlJ. This is 



because the effect of Pauli exclusion principle cancels out (in analogy with Boltzmann 
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picture with elastic scattering). Substituting p(E) = g(E)f(E — E F ), where f(E — E F ) 
is the Fermi occupation function, one obtains 

H = \g{E F ) VF (A.6) 

where v F should be calculated at the Fermi energy. 

Finally, the microcanonical version of the FD relation is 

- = lwm lg{E} ^ (A - 7) 

The subscript E indicates that both is E and /i E are evaluated locally around some 
energy E. 



Appendix B. Cross correlations I 

In this Appendix we introduce two proofs of Eq. (p4|) . The first is a formal argument, 
while the second is a more physically appealing argument. The formal argument is 
as follows: Eq.(p"2|) is an exact result which can be written using obvious abstract 
matrix notation as v = D-f E D. Let D = D\ + L>2- If D2 is a special deformation 
then by definition -D27E-D2 = 0. But this can be true only if D2 belongs to the 
kernel (nullspace) of the matrix 7 E) hence we have 7 E Z?2 = 0. Therefore we have also 
Z)i7 E L*2 = for any D\, which is precisely the statement of Eq. (p4[) . 

Now we present the alternative physically appealing argument. Consider two 
noisy signals Tit) and Q{t). We assume that (Tit)) = (Q(t)) = 0. The angular 
brackets stand for an average over realizations. The auto-correlations of Tit) and Q(t) 
are described by functions C F (r) and C g (t) respectively. We assume that both auto- 
correlation functions are short-range, meaning no power-law tails (this corresponds to 
the hard chaos assumption of this paper), and that they are negligible beyond a time 
r„. We call a signal 'special' if the algebraic area under its auto-correlation is zero. 
The cross-correlation function is defined as 

c p ,g(t) = (r(t')g(t")), T = t'-t". (b.i) 

We assume stationary processes so that the cross-correlation function depends only 
on the time difference r. We also symmetrize this function if it does not have r 1— » — r 
symmetry. We assume that C F _ G (r) is short-range, meaning that it becomes negligibly 
small for |r| > r c . We would like to prove that if either J-(t) or Q(t) is special then 
the algebraic area under the cross-correlation function equals zero. 

Consider the case where Tit) is general while Q(t) is special. The integral of 
C f (t) will be denoted by v. Define the processes 

X(t) = f T{t')dt' (B.2) 
Jo 

Y(t) = f Git")dt". (B.3) 
Jo 

From our assumptions it follows, disregarding a transient, that for t r c we have 
diffusive growth (X(£) 2 ) w vt. However since Yit) is a stationary process |32|, 
(F(i) 2 ) const. Therefore for a typical realization we have < const x V vt 

and < const. Consequently, without making any claims on the independence of 
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X(t) and Y(t), we get that (X(t)Y(t)) cannot grow faster than const x \fvi. Using 



the definitions (JB.2|), (B.3) and (B.l) we can write 



Cjo(T)dT= W)y( f )) w co^t^ (R4) 

> t y/t 

where the limit t — > oo is taken. Thus we have proved our assertion. 
Appendix C. Cross correlations II 

In this section we further discuss some features of the cross-correlation function. For 
the purpose of presentation we we would like to view the time as an integer variable 
t = 1,2, 3.... One may think of each instant of time as corresponding to a bounce. 

Let us assume that we have functions /(s) and g(s), and a time-sequence 
(si, S2, S3, ...). This gives two stochastic-like processes (JFi, JF 2 . JF 3 , ...) and 
(Gi, G2, G3, ••■)■ The cross correlation of these two processes is defined as follows: 

^(i-j) = W) = (/Ws('i)) (CI) 
It is implicit in this definition that we assume that the processes are stationary, so the 
result depends only on the difference r = {i — j). The angular brackets stand for an 
average over realizations of s-sequences. 

If the sequences are ergodic on the s domain, then it follows that 



(F> = Jf(s)ds 

(G) = Jg{s)ds 

a, G (0)= [f(s)g(s)ds (C.2) 



The t / cross-correlations requires information beyond mere ergodicity. In case 
that the s sequence is completely uncorrelated in time we can factorized the averaging 
and we get C FjG (r ^ 0) = (T) x (£). If (T) = then 

C f , g (t^0)=0 (C.3) 

irrespective of (Q) . 

However, we would like to define circumstances in which Eq.( |C.3| ) is valid, even 
if the s sequence is not uncorrelated. In such case either the T or the Q may possess 
time correlations. (Such is the case if Q is 'special'). So let us consider the case where 
the T sequence looks random, while assuming nothing about the Q sequence. By the 
phrase 'looks random' we mean that the conditional probability satisfies 

Prob(^|sj) = Prob(^) for any i ^ j (C.4) 



Eq. (C.3) straightforwardly follows provid ed (! F) = 0, irrespective of the g(s) involved. 



Given /(s), the goodness of ass umption ( |C.4| ) can be actually tested. However, it is 



not convenient to consider (C.4) as a practical definition of a 'normal' deformation. 
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